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We review recent progress in the development of an approach valid to any order which unifies the fixed-order 
DGLAP evolution of fragmentation functions at large x with soft-gluon logarithmic resummation at small x. 
At leading order, this approach, implemented with the Double Logarithmic Approximation, reproduces exactly 
the Modified Leading Logarithm Approximation, but it is more complete due to the degrees of freedom given 
to the quark sector and the inclusion of the fixed-order terms. We find that data from the largest x values to 
the peak region can be better fitted than with other approaches. In addition, we develop a treatment of hadron 
mass effects that leads to additional improvements at small x. 



1. Introduction 

The perturbative approach to Quantum Chromody- 
namics (QCD) is believed to solve all problems within 
its own limitations provided the correct choices of the 
expansion variable and the variable(s) to be fixed are 
used. However, perturbative QCD (pQCD) currently 
has the status of being a large collection of seem- 
ingly independent approaches, since a single unified 
approach valid for all processes is not known. This is 
problematic when one wants to use a range, qualita- 
tively speaking, of different processes to constrain the 
same parameters, e.g. in global fits. What is needed 
is a single formalism valid over the union of all ranges 
that the various pQCD approaches allow. This uni- 
fication must be consistent, i.e. it must agree with 
each approach in the set, when the expansion of that 
approach is used up to the order being considered. 

In this presentation, we review recent progress in 
establishing a generalized formalism for the evolution 
of fragmentation functions (FFs) from the smallest 
x values probed so far way up to unity Q. After 
explaining how to treat soft-gluon and hadron-mass 
effects in Sees. |2 and |3 respectively, we study their 
phenomenological implications in Sec.0| Our conclu- 
sions are contained in Sec. |S] 



2. Resummation of soft-gluon 
logarithms 

The current optimum description of single-hadron 
inclusive production is provided by the QCD par- 
ton model, which requires knowledge of the FFs 
D^(x,Q 2 ), each corresponding at leading order (LO) 
to the probability for the parton a produced at short 
distance 1/Q to form a jet that includes the hadron h 
carrying a fraction x of the longitudinal momentum of 
a. From now on we omit the label h and put the FFs 
for all a into the vector D(x,Q 2 ). The evolution of 
the FFs in the factorization scale Q 2 is also required, 
which at large and intermediate x is well described 



by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equation |3J, given at LO by 



d\nQ 



;D{x,Q 2 ) 



dz 



z a s (Q 2 )P^(z)D v 



Q~ 



(i) 

where P(°'(z) are the LO splitting functions cal- 
culated from fixed-order (FO) pQCD. We define 
a s = a s /(2ir), which at LO obeys a s (Q 2 ) = 
l/[/? ln(Q 2 /A 2 QCD )], where (3 = (11/6)0* - 
(2/3)Tr71/ is the first coefficient of the beta function 
and Aqcd is the asymptotic scale parameter of QCD. 
For the color gauge group SU(3), the color factors ap- 
pearing in this presentation are Cf = 4/3, Ca = 3, 
and Tji = 1/2; nf is the number of active quark fla- 
vors. 

On the other hand, at small x, the FO approxi- 
mation fails due to unresummed soft-gluon logarithms 
(SGLs). The most singular SGLs are the double loga- 
rithms (DLs), which give a z — > singularity in the LO 
splitting function a s P^(z) of the form a s /z. These 
DLs occur at all orders in the FO splitting function, 
being generally of the form (l/z) (a s In z) 2 (a s In 2 z) r 
for r = — l,...,oo, implying that, as x decreases, 

Eq. (TJ3) becomes a poor approximation once ln(l/a;) = 
r /2 

0(a s ). An improvement of the small- X description 
should be obtained by accounting for all DLs, which is 
provided by the Double Logarithmic Approximation 
(DLA)fia, 



d\nQ 



where A 
while 



:D(x,Q 2 ) 



1 dz 2C A 



Az dl =« 2 

z z 



x D(-,Q' 

. z 



a s (Q 2 



(2) 



for valence-quark or non-singlet FFs, 



A = 



n 2< ^ f 

U C A 

I 



(3) 



when D = (D^.D g ), where D E = (l/n f ) EgU( D Q + 
Dq) is the singlet FF. The Modified Leading Loga- 
rithm Approximation (MLLA) 0, 0, improves the 
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description over the DLA only by including a part of 
the FO contribution that is known to be important 
at small x. With certain qualifications g|, the MLLA 
leads to a good description of all data down to the 
smallest x values. 

However, what has been lacking is a single approach 
which can describe data from the largest to the small- 
est x values. Such an approach can be simply and 
consistently constructed from the knowledge provided 
by the DLA and the DGLAP equation. The result re- 
produces the FO result when expanded in a s , while all 
small-w singularities are resummed into a non-singular 
expression in Mellin space (see later) . As shown in de- 
tail in Ref. this approach is unique up to higher- 
order terms in the FO series and up to less singular 
SGLs, e.g. the single logarithms (SLs). Such a formal- 
ism may be constructed simply by evolving according 
to Eq. Q , but with the splitting function modified as 



s P^(z)^P Bh (z,a s ) + . 



(4) 



where P DL (z,a s ) contains the complete resummed 
contribution to the splitting function to all orders from 

the DLs, while a s P^°\z) is the remaining FO contri- 
bution at LO. It is obtained by subtracting the LO 
DLs, already accounted for in P DL , from a s P(°\z) to 
prevent double counting. We work to LO, since the 
less singular SGLs that occur at next-to-leading order 
(NLO) are not known to all orders. However, we leave 
the SL at LO unresummed, since it is not singular as 
z — ► 0. We now use Eq. J3J) to gain some understand- 
ing of P DL . For this, we transform to Mellin space, 
defined by 



f(uj)= / dxx u f(x). 



(5) 



Upon Mellin transformation, Eq. © becomes 



ui + 2 



d\nQ 2 J dlnQ 2 
2C A a s (Q 2 )AD{Lj,Q 2 ) 



D(u,Q 2 ) 



(6) 



Substituting Eq. |@} with the FO term a s P (z) ne- 
glected in Eq. taking the Mellin transform, 

D(uj,Q 2 ) = P BL (u;,a s (Q 2 ))D(u J ,Q 2 ), (7) 



d\nQ 2 

and inserting this in Eq. JSJl, we obtain 

2(P DL ) 2 + uP DL - 2C A a s A = 0. (8) 
We choose the solution 

P dl (uj, a.) = j (-co + ^u 2 + 16C A a s ) , (9) 

since its expansion in a s yields at LO the result 

a s P^\co,a s )=( Q (10) 



P 



DL 



which agrees with the LO DLs from the literature 0] • 
Equation (fjljl contains all terms in the splitting func- 
tion of the form (a s /uj)(a s /uj 2 ) r+1 , being the DLs in 
Mellin space, and agrees with the results of Refs. @,@] . 
We now return to x space, where Eq. (Fj|J) reads 

(z,a s ) = A f^ Ji (aVCaVs^-) , (11) 
zln(l/z) \ z J 

with Ji being the Bessel function of the first kind. 

To summarize our approach, we evolve the FFs 
according to Eq. Q, but with the replacement of 
Eq. (@J, where P DL (z,a s ) is given by Eq. (fTTf) . and 

a s P^°\z) is given by a s P^(z) after the terms pro- 
portional to a s /z have been subtracted. 

To extend NLO calculations to small x p , the com- 
plete resummed DL contribution given by Eq. (|llf) 
must be added to the NLO splitting functions. These 
contain SGLs belonging to the classes m = 1, ... ,4, 
which must be subtracted. Note that the NLO m = 1 
term is accounted for by the resummed DL contribu- 
tion. The m — A term is a type p = term, and hence 
does not need to be subtracted. 

Before we outline the phenomenological investiga- 
tion of our approach, we note that it is more complete 
than the MLLA, which can be shown as follows. With 

a s P^ (z) accounted for, Eq. J()J) is modified to 



dlnQ 2 j dlnQ 2 
2C A a s {Q 2 )AD{oj,Q 2 ) 

x p ( V)r,( w ,Q 2 )l , 



D(lo,Q 2 ) 



rflnQ 2 



a s (Q 2 ) 
(12) 



up to terms which are being neglected in this presen- 
tation and which arc neglected in the MLLA. If we 

approximate a s P^°\uj) by its SLs, defined at LO to 
be the coefficients of uj°, equal to those in a s P^(u)), 



p sl (°)H = ( 



-3C F 



(13) 





\T R n t -±C A - jT R n f 

and apply the approximate result that follows from 
the DLA at large Q, 



D - = 9Ld 

q q ~ C A 9 ' 



(14) 



which can be derived, e.g., from Eq. Q11JI. the gluon 
component of Eq. Ijl2|l becomes precisely the MLLA 
differential equation. Therefore, we conclude that, 
since we do not use these two approximations, our 
approach is more complete and accurate than the 
MLLA. 



3. Incorporation of hadron mass effects 

We now incorporate hadron mass effects into our 
calculations, using a specific choice of scaling variable. 
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For this purpose, it is helpful to work with light cone 
coordinates, in which any four-vector V is written in 
the form V = (V + ,V-,V T ) with V ± = (V°±V 3 )/V2 
and V T = (V 1 , V 2 ). In the center-of-mass (CM) 
frame, the momentum of the electroweak boson takes 
the form 



y/2' V2 



,0 



(15) 



In the absence of hadron mass, x p (whose definition 
x p = 2p/y/s applies only in the CM frame) is identical 
to the light cone scaling variable n = p^ /q + . However, 
the definition x v — 2p/s/s applies only in the CM 
frame. So, 77 is a more convenient scaling variable for 
studying hadron mass effects, since it is invariant with 
respect to boosts along the direction of the hadron's 
spatial momentum. Taking this direction to be the 
three-axis, and introducing a mass rrih for the hadron, 
the momentum of the hadron in the CM frame reads 



Ph 



= ( vVl m l 







(16) 



Therefore, the relation between the two scaling vari- 
ables in the presence of hadron mass is 



mi 

= Tl[l--± 

ST] 



(17) 



Note that these two variables are approximately equal 
when <C x p ^/s, i.e. hadron mass effects cannot be 
neglected when x p (or 77) are too small. 

In the leading-twist component of the cross section 
after factorization, the hadron h is produced by frag- 
mentation from a real, massless parton of momentum 



k = 



2L 
y 



,0,0 



(18) 



The + component of everything other than this parton 
and of everything produced by the parton other than 
the observed hadron h must be positive, implying that 
y > 77 and y < 1, respectively. As a generalization of 
the massless case, we assume the cross section we have 
been calculating is (da / 'drj)(rj, s), i.e. 

= [ l ±*to,e,(?)DK(A, (19) 
dv J v y dy \y ) 

which is related to the measured observable 
(da / dx p )(xp, s) via 



da 



1 



da 



dx p ^ p ' ^ 1 + m 2 i /[srj 2 (x p )] drj 



(V(x p ),s). (20) 



Note that the effect of hadron mass is to reduce the 
size of the cross section at small x v (or n). 

The data we are studying are described by the sum 
of the production cross sections for each light charged 



hadron species, being the charged pion, the charged 
kaon, and the (anti)proton, whose masses (140, 494, 
and 938 MeV, respectively) are substantially different. 
Therefore, separate FFs and hadron masses for each 
of the three species are needed. Clearly, so many free 
parameters would not be constrained by these data. 
However, since most of the produced particles are pi- 
ons, it is reasonable to take the final-state hadrons to 
have the same mass, so that in our approach it is rea- 
sonable to use a single hadron mass parameter m/j, 
whose fitted result should be closer to the pion mass 
than the proton mass, i.e. around 300 MeV. Any sig- 
nificant deviation from this value would imply that 
some physics has not been accounted for. 



4. Numerical analysis 

In a first test of our approach, we compare its effects 
on fits of quark and gluon FFs to e + e _ data with the 
standard FO DGLAP evolution. We use normalized 
diffcrential-cross-section data for the process e + e~ — * 
(7, Z) — > h + X, where h is a light charged hadron 
and X is anything else, from TASSO at y/s — 14, 
35,44 @, and 22 GeV [13 , MARK II [HI and TPC 
Hi at 29 GeV, TOPAZ at 58 GeV [H, ALEPH_E1, 
DELPHI L3 [ll, OPAL [Hf, and SLC [Hjat 
91 GeV, ALEPH [lland OPAL ^ at 133 GeV, 
DELPHI at 161 GeVJ2l|, and OPAL at 172, 183, 189 
0, and 202 GeV |23f~~We place a small-i cut 
our data of 



^ = l n I<ln^ ; 
x 2M 



(21) 



where M is a mass scale of 0(1 GeV). We fit the 
gluon FF, D g (x, Qq), as well as the quark FFs, 

D dsb {x,Ql) - ^[D d {x,Ql)+D s {x,Ql) 

+ D b (x,Q 2 )] , 
D uc {x,Ql) = l - [D u {x,Ql) + D c {x,Ql)} , (22) 



where Qo — 14 GeV. Since the hadron charge is 
summed over, we set Dq = D q . For each of these 
three FFs, we choose the parameterization 

D(x,Ql) = A^exp(-cln 2 x)a; a (l-x) /3 , (23) 

which, at small x, is a Gaussian in £ for c > with 
center positive in £ for a < as is found to be the 
case, while it reproduces the standard parameteriza- 
tion (i.e. that without the exp(— cln 2 a;) factor) used 
in global fits at intermediate and large x. We use 
Eq. (|14f> to motivate the simplification 



Cdsb 



'91 



OLdsb 



(24) 
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to our parameterization. We also fit Aqcd, so that 
there is a total of 9 free parameters. Since we only 
use data for which y/s > mb, where TO& « 5 GeV is 
the mass of the bottom quark, and since Qo > m^, we 
take rif — 5. While the precise choice for n/ does not 
matter in the DLA, calculations in the FO approach 
strongly depend on it. 

Since the theoretical error in our LO approach is 
larger than that at NLO, some of our results should 
be interpreted somewhat qualitatively. In particular, 
the result for Aqcd has a theoretical error of a factor 
ofO(l). 

Without the (1 — x)@ factors, Eq. lfl"4")) implies that 



N,, 



N, 



dsb 



Ca 



(25) 



Since the (1 — x)@ factors are important at large x, 
where Eq. (|14(l is no longer valid, we do not impose 
Eq. J23J, but rather test its validity over the whole x 
range. 



Table I Parameter values for the FFs at Qo = 14 GeV 
parameterized as in Eq. 12311 from a fit, with Xdf = 3-0, 
to all data listed in the text using DGLAP evolution in 
the FO approach to LO. The fit also yields 
Aqcd = 388 MeV. 



FF N 


P 


a 


c 


g 0.22 


-0.43 


-2.38 


0.25 


u + c 0.49 


2.30 


[-2.38] 


[0.25] 


d+s + b 0.37 


1.49 


[-2.38] 


[0.25] 




Figure 1: Fit to data as described in Table Some of 
the data sets used for the fit are shown, together with 
their theoretical predictions from the results of the fit. 
Data to the right of the vertical dotted lines are not 
used. Each curve is shifted up by 0.8 for clarity. 



We first perform a fit to all data sets listed above us- 
ing standard LO DGLAP evolution, i.e. Eq. (JTJ with- 
out the replacement in Eq. (0J. We fit to those data 
for which Eq. J2TJ) is obeyed with M = 0.5 GeV. This 
gives a total of 425 data points out of the available 492. 
We obtain Xbf ~ ^.0 (or 2.1 after subtraction of the 
contribution to x 2 from the TOPAZ data, which is the 
only data set whose individual Xbf ^ s greater than 6), 
and the results are shown in Table [I] and Fig. It is 
clear that FO DGLAP evolution fails in the descrip- 
tion of the peak region and shows a different trend 
outside the fit range. The exp(— cln 2 :r) factor does 
at least allow for the fit range to be extended to x 
values below that of x = 0.1, the lower limit of most 
global fits, to around x = 0.05 (£ = 3) for data at the 
larger s/s values. Note that (3 g is negative, while kine- 
matics require it to be positive. However, this clearly 
does not make any noticeable difference to the cross 
section. 



Table II Parameter values for the FFs at Qo = 14 GeV 
parameterized as in Eq. I23H from a fit, with Xdf = 2.1, 
to all data listed in the text using DGLAP evolution in 
the DLA-improved approach to LO. The fit also yields 
Aqcd = 801 MeV. 



FF 



N 



g 1.60 5.01 -2.63 0.35 
u + c 0.39 1.46 [-2.63] [0.35] 
d + s + b 0.34 1.49 [-2.63] [0.35] 



20 



15 - 



• OPAL 202 
« OPAL 1 72 

a DELPHI 161 

• OPAL 133 
. OPAL 91 

□ TASSO 44 
« TASSO 35 

• TPC 29 

■ TASSO 1 4 



S 10 - 



5 - 



J0 





1 2 3 4 5 6 7 

\ 

Figure 2: Fit to data as described in Table ITT1 

Now we perform the same fit again, but using our 
approach, i. e. Eq. with the replacement in Eq. Q , 
for the evolution. The results are shown in Table ITTl 
and Fig. |2 We obtain Xdf = 2.1 (or 1.4 without 
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the TOPAZ data, the individual Xdf f° r eacn remain- 
ing data set being less than 3), a significant improve- 
ment relative to the fit above with FO DGLAP evo- 
lution. This should also be compared with the fit to 
the same data in Ref. (24|, where DL resummation 
is used within the MLLA but with neither FO terms 
nor quark freedom {i.e. Eq. (|f 4|l is imposed over the 
whole x range) and Xdf = 4-0 is obtained. The data 
around the peak is now much better described. The 
energy dependence is well reproduced up to the largest 
s/s value, 202 GeV. At y/s = 14 GeV, the low-x de- 
scription of the data is extended from x = 0.1 in the 
unresummed case down to 0.06 in the resummed case, 
and from x = 0.05 to 0.005 at yfi = 202 GeV. 

We conclude that, relative to the MLLA, the FO 
contributions in the evolution, together with free- 
dom from the constraint of Eq. (|14f> . make a sig- 
nificant improvement to the description of the data 
for £ from zero to just beyond the peak. The value 
Aqcd ~ 800 MeV is somewhat larger than the value 
480 MeV, which we obtain from a DGLAP fit in the 
large- a; range (x > 0.1). Had we made the usual 
DLA (MLLA) choice Q = y/s/2 instead of our choice 
Q = ^/s as is done in analyses using the DGLAP equa- 
tion, we would have obtained half this value for Aqcd- 
In an analysis at NLO, the choice of Q is less relevant, 
since the theoretical error on Aqcd is smaller. In ad- 
dition, while the DL resummation greatly improves 
the description around the peak, it is still not per- 
fect. N g exceeds the value expected from Eq. ifHEjl by 
a factor of about 2. However, note that N g is weakly 
constrained, since the gluon FF couples to the data 
only through the evolution, requiring gluon data to 
be properly constrained. These problems are related 
to the worsening of the description of the data on mov- 
ing beyond the peak, since fits in which the cuts are 
moved to larger £ values give an increase in Aqcd and 
N g , as well as in Xdf- 

Figure[21is repeated in Fig. [3] to show more clearly 
the good quality of the fit at intermediate and large x. 
A couple of points at x = 0.9 are not well described, 
although the data here are scarce and have larger er- 
rors. This shows that the DL resummation does not 
deteriorate the quality of the FO description at large 
x. This is further confirmed by noting that, for the 
above fit in our approach, the subsample of the data 
at x > 0.1 yields Xdf = 3-0, which is close to the 
result Xdf = 2.9 obtained from fitting to the same 
large-x data using only the FO approach. 

We now perform the last fit again, but with m,h in- 
cluded in the list of free parameters. We obtain the 
results in Table UTTI and Fig.^J and the values Aqcd = 
399 MeV and m h = 252 MeV. Treatment of hadron 
mass effects renders the value of Aqcd obtained in the 
fit with DL resummation more reasonable and slightly 
improves the latter, yielding Xdf = 2-0- We conclude, 
therefore, that to improve the large-£ description and 
to achieve a reasonable value for Aqcd, both DL re- 



. OPAL 202 
a OPAL 172 
s DELPHI 161 

■ OPAL 133 

■ OPAL 91 

□ TASSO 44 
« TASSO 35 
• TPC 29 
. TASSO 14 




Figure 3: As in Fig. [5] but with the cross section on a 
logarithmic scale versus x. Each curve, apart from the 
lowest one, has been rescaled relative to the one 
immediately below it by a factor of 5 for clarity. 



Table III Parameter values for the FFs at Qo = 14 GeV 
parameterized as in Eq. II23II from a fit, with Xdf = 2.0, 
to all data listed in the text using DGLAP evolution in 
the DLA-improved approach to LO incorporating hadron 
mass effects. The fit also yields Aqcd = 399 MeV and 
m h = 252 MeV. 



FF 



N (3 



g 1.59 7.80 -2.65 0.33 
u + c 0.62 1.43 [-2.65] [0.33] 
d + s + b 0.74 1.60 [-2.65] [0.33] 




Figure 4: Fit to data as described in Table ITTT1 



summation and treatment of hadron mass effects are 
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required. 

5. Conclusions 

In conclusion, we have proposed a single unified 
scheme which can describe a larger range in x than 
either FO DGLAP evolution or the DLA, by imple- 
menting SGL resummation in the standard DGLAP 
formalism in a minimal and unambiguous way through 
any order of perturbation theory. Our present limita- 
tion to LO is caused by the lack of knowledge on SGL 
resummation necessary to match NLO DGLAP evolu- 
tion. The actual inclusion of higher orders, although 
conceptually straightforward in our scheme, requires 
a separate analysis, which lies beyond the scope of 
Ref. [l| and is left for future work. Further improve- 
ment in the small- x region can be expected from the 
inclusion of resummed SLs. Our scheme allows for a 
determination of quark and gluon FFs over a wider 
range of data than previously achieved, and should be 
incorporated into global fits of FFs such as that in 
Ref. [25|, since the current range of 0.1 < x < 1 is 
very limited. 

Since FFs and their DGLAP evolution are univer- 
sal, our approach should be expected to also improve 
the description of inclusive hadron production from 
reactions other than e + e~ annihilation, such as those 
involving at least one proton in the initial state. 
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